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We present an analytic theory for the energy spectrum of 
a two-leg spin ladder doped with two holes. Starting from a 
pseudo-fermion-bond-boson representation of the correspond- 
ing ti,2 — Ji,2 Hamiltonian we apply a diagrammatic approach 
adapted to the limit of strong rung coupling, which includes 
both, the coupling of holes to the spin background as well 
as the two-hole interactions. The two-hole spectrum is cal- 
culated and the formation of bound states is discussed. Ad- 
ditionally the evolution of the spin gap of the ladder upon 
doping is analyzed. A comparison with existing exact diago- 
nalization data is presented and good agreement is found. 



I. INTRODUCTION 

Low dimensional quantum-spin systems with local mo- 
ments arranged in ladder-type geometries have attracted 
considerable interest in recent years. From a theoretical 
point of view spin ladders are of interest as. a bridging 
step from one- to two-dimensional systemsEJ. From an 
experimental Boint of vieW|*"arious spin ladder materials, 
e.g. SrCu2 03l3 or CaV2 05H, have been analyzed, allow- 
ing to gauge theoretical findings against real systems. 

Spin-1/2 two-leg ladders display a quantum disor- 
dered ground state with exponentially decaying spin- 
correlations and a gapful magnetic excitation spectrum. 
Apart from magnetic excitations the dynamics of charge 
carriers doped into the quantum disordered ground state 
is of interest, in particular because of the discovery 
of superconductivitycl under pressure in the spin-ladder 
compound Sri4_2;Ca2;Cu2404i3. The properties of a sin- 
gle hole doped into the i^l^ ladder have been the-subp 
ject of various numcricalQ^li^ and analytical studiesliBilJ. 
Moreover, in order to understand two-hole interactions 
as well as pairing correlations several numerical investi- 
gations have been performed focusinffi-pp-|Various topics 
like the fopaation of bound hole-pairsRljI^l, exojmd state 
propertiesE£tEa-ami excitation spectrumBtrEJiH3, coi; 
tion functional3'0 and superconducting fluctuationsc 
Regarding analytical approaches however, only a re- 
stricted set of pxesults is available like hifih, order se- 
ries expansion£3, meaaj-field approachesE3c3, effective 
Hamiltonian mappingsEy, perturbative renormaliaation 
group analysis^, exact solutiop-pf a related modelE2l and 
recurrent variational approachcJ. In particular, the two- 
hole spectrum remains an interesting issue with many 
open questions. 



In this work we present an analytic theory of the energy 
spectrum of a two-leg spin ladder doped with two holes. 
Our approach allows for a direct understanding of the 
relevant processes. Moreover, while showing|-r( 
good agreement with earlier numerical studiesO 
on small systems, our method may equally well be ap- 
plied to large systems close to the thermodynamic limit. 

For the two-leg spin ladder we consider a i — J-model 
with hopping- and exchange-integrals along the rungs 
and chains of fig. |l|, i.e. ti, Ji and t2,J2 respectively. 
The Hamiltonian reads 
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where Si.„ are spin-1/2 operators on site i of rung n and 

^i,n,a = [ci,n,<T(l " "i.ri.-a)]^^^ are fcrmion-opcrators for 
spin a and site i, n projected onto the space of no double 
occupancy. 

The paper is organized as follows. In the next sec- 
tion we briefly restate the bond-operator description of 
the spin dynamics of the undoped ladder. In section III 
a mapping of the t-J Hamiltonian in the two-hole sec- 
tor onto an interacting fermion-boson-model is described. 
Section pM presents a diagrammatic analysis for the en- 
ergy spectrum. In section M results are presented and 
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FIG. 1. Two-leg t-J-ladder. 
denote the sites on the rung. 



n+l 
n labels the rungs. 
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compared to existing numerical data. Finally conclusions 
will be given. 



II. SPIN DYNAMICS OF THE UNDOPED 
LADDER 

To describe the spin excitations of the two-leg lad- 
der in the sector of no holes jaicj_use the well established 
bond-operator representationlUa of dimerized spin-1/2 
systems. Here we briefly recapitulate this method. The 
eigenstates of the total spin on a single rung occupied by 
two electrons are a singlet and three triplets. These can 
be created by the bosonic bond operators sj^ and tj, with 
a = x,y, z acting on a vacuum |0) by 
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where the first (second) entry in the kets refers to site 1(2) 
of the rung n of Fig. n^. On each site we have [s,s^] — 1, 



The bosonic Hilbert 



[s^'^4''] =0, and [to ,tl] =Sc,f3. 

space has to be restricted to either one singlet or one 

triplet per site by the constraint 
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Expressing the spin part Hj of the Hamiltonian (ffl) 
by the bond operators yields an bose gas of singlets and 
triplets with two-particle interactions mediated by the 
inter- rung coupling J2. At J2 = these interactions 
vanish leaving a sum of purely local rung Hamiltoni- 
ans, which lead to a product ground-state of singlets 
localized on the rungs and a set of 3''^-fold degenerate 
triplets. For finite J2 the inter-rung interactions can be 
treated approxiiaatelv-by a linearized Holstein-Primakoff 
(LHP) approachtil'E3T2a. The LHP method retains spin- 
rotational invariance and reduces Hj to a set of three 
degenerate massive magnons 
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The '[] '-bracketed sign on the rhs. in (H) refers to the 
quantity v on the Ihs.. The spin-gap A = minjwfc} re- 
sides at fc = TT with A = y/Ji ~ 1J\J2- Note that be- 
cause of (H) the ground state ID), which is defined by 
7a,fc|D) = 0, contains quantum-fluctuations beyond the 
pure singlet product-state. To leading order the [disper- 
sion LOk is identical to perturbative expansionsElo. Be- 
yond the LHP approach triplet interactions and more 
elaborate consideration of the constraint (||) lead to a 
renormalizatiojjLofjiiife and the formation of multi-magnon 
bound statescirc3il3. However, for J2 <C J\ these renor- 
malizations can be neglected. 



III. TWO-HOLE HAMILTONIAN 

Li the two hole sector a mapping of the Hamiltonian (Q) 
onto a rung basis requires for operators which describe 
singly occupied as well as empty rungs in addition to 
the singlet and triplet operators (g) of the undoped spin 
system. We label a rung occupied by aj-sJngle-hole by 
introducing an additional pseudo-fermionEj (holon). Le. 
instead of the rung being in one of the states given by 
(0) it can also be in the states 



<«.J0) = \30)r 



(9) 



where the l.h.s. denotes the vacuum, i.e. |0), with a 
single rung at site n in a one-hole state of spin a with 
i = 1,2 referring to the two positions available for the 
hole on the bond. Moreover two holes on the same rung 
are created by an additional boson 
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where now the l.h.s. refers to the vacuum |0) with one 
hole-pair at site n. Both, the single-hole operators a^ ^ cr 
and the hole-pair operators &„ are required to obey the 
standard fermionic and bosonic operator algebras. 

In the two-hole sector each rung necessarily is in one of 
the states (||), (g) or (nO|). As a consequence an extended 
hard-core constraint (^ has to be fulfilled 

with a summation over repeated spin-indices implied. 

Creation of a physical hole in the (half-filled) ground 
state \D) of the spin-system is described byc3 

Cj",n,cr = —7=[Cl-j^_(PjPa-Sn <*z,n) 

+<^„JPo^tx,n + ity,n) ] + ^\ai,«,<T (12) 

where pj = +(— ), j = 2(1) for j ~ 1(2) and Pa = +(— ), 
W ^l (t) for a —^ (l). In terms of ( [l^ ) the creation of a 
physical hole can be interpreted as the removal of a spin- 
rung followed by the creation of a single-hole state or the 



removal of a single-hole state followed by the creation of 
a hole-pair state. The particular linear combination of 
the a}-^^ s„, and t^„^ operators on the r.h.s. ensures 
that the total spin is S* = 1/2 and 5*^ = ±1/2. Using the 
constraint (O) and the bosonic (fermionic) commutation 
relations for the &„ (Qj.n.o-) it is straightforward to show, 
that the r.h.s. of ( [l2|) indeed satisfies the usual Hubbard- 
operator algebra. The pseudo-particle representation of 
the spin operators reads 
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where t"^, are Pauli-matrices. 



Inserting (y_2|) and (|1^) into Hamiltonian (|1|) various 
processes are found: First, single- hole processes like (i) 
intra-rung hopping, (ii) inter-rung hopping, and (iii) ex- 
change scattering of the pseudo-fermions. Second, two- 
hole processes result like the (iv) decay of a hole-pair on 
a rung and a singlet (triplet) into two separated pseudo- 
fermions, (v) exchange scattering between two pseudo- 
fermions on neighboring rungs and (vi) an on-site bind- 
ing energy of the hole-pair, due to the additional singlet 
in the system as compared to a situation with two holes 
occupying separate rungs. Finally, three-hole processes 
occur which involve one hole-pair and one pseudo-fermion 
creation operator. These can be discarded in the two-hole 
sector. Processes (i) - (iii) constitute the Hamiltonian Hj 
of the one-hole sector which has been discussed in ref.EJ. 
For the sake of completeness Hi is listed in appendix 
The remaining processes (v) and (vi) are contained 
in Hh while Ht represents process (iv) of the the final 
Hamiltonian which reads 







H = Hi + Ht + Hj 
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where j ^ 2(1) for j = 1(2) and Pa = +(-), '^ ==i (T) 

for Cr =t (i). S5f,_„ = hT.a,,a2 4,m,a^'^<yl,'r2aj^n^^., and 

n° is the number operator of fermion aj^„. Following 
the LHP approximation we have replaced the singlet op- 
erators by C-numbers, i.e. s. Within the LHP approach 



s = 1 to lowest order. Regarding the charge dynam- 
ics we improve upon this by requiring that the conden- 
sate density of the singlet is determined by satisfying the 
hard-core constraint (^ on the average, i.e. by setting 
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In principle this equation can be used to determine the 
LHP magnon-dispersion selfconsistently which however 
we will refrain from here. 

After Fourier transforming (HJlq) and inserting the 
Bogoliubov representation (|^) we get 
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[cos {q - q') +COS {k-q- q')] (23) 



and the (anti)bonding linear combinations of single-hole 

states 
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which have momentum 0(7r) i.e. even (odd) parity in 
rung direction, where the -l-(— ) belongs to 0(7r). In ([2l|), 
and as we focus on the limit of J2 «C Ji , where the LHP 
provides for a controlled approximation of the spin sys- 
tem, we have dropped the exchange scattering ex J2 from 
Hbof{^. 

In the two-hole sector the hard-core constraint (|ri|) is 
of particular importance in order to suppress double oc- 
cupancies of a single rung by two a- type pseudo-fermions. 
This cannot be achieved by enforcing ( |l l| ) only on the av- 
erage, rather an auxiliary intra-rung hard-core potential 
has to be added to the Hamiltonian ([lq-pi|) 
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with C/q —^ cxD. A similar hard-core potential, howeveLto 
study the two-triplet sector has been used in refs.Cjifj. 



IV. TWO-HOLE SPECTRUM 

The primary goal of this work is to evaluate the mo- 
mentum dependent energy spectrum of a two-hole state 
on the ladder. While the physical single-hole excitations 
are composite objects, c.f. (fsh, their deviations from 
the pseudo fermions are small for J2 <C Ji, i.e., for a 
small triplet |_dcnsity induced by quantum fluctuations 
(see also rcf.Efl). Therefore, as a first step towards our 
goal we approximate the two-hole energy spectrum by 
that of two pseudo-fermions based on the limit of strong 
intra-rung exchange. This leaves (i) a two-particle prob- 
lem to be solved, with however (ii) strongly renormalized 
one-particle excitations. On the bare two-particle level 
(i) can be achieved exactly by summing the T-matrix 
incorporating the interactions Hs, Ht and U. For an 
approximate account of (ii), we will sum the T-matrix 
using renormalized onc-particlc pseudo-fcrmion Green's 
functions 

Gfe„.(fc,i) = -^Qit){D\{au^^uAtWk,.u,<r}\D)■ (26) 

The main source of renormalization for a single pseudo- 
fermion is known to be multi-triplet emission which can 
be accounted for by a selfconsistent resummatien of all 
noncrossing n-loop graphs to the self energyllJ. This 
is equivalent to the selfconsistent-Born-approximation 
(SCBA) which has been emplpyed for the one- and 
two-dimensional t-J-modelcZTE^. On the ladder the 
SCBA leaves a finite quasi-particle residue to the pseudo- 
fermion propagator, which is a consequence of the spin 
gap. Yet, the spectra acquire a substantial incoherent 
part as well as renormalization of the quasi-particle dis- 
persion (for details see ref.Eil). 

To evaluate the T-matrix we note that the two-leg 
ladder is symmetric with rcspcct_to reflections at a 
plane perpendicular to the rungSLj. This leaves the 
corresponding parity a good quantum number. The 
the (anti)bonding single-hole states aL^oncrl^) ^"^^ *^^ 
even(odd) parity, the hole-pair state b'l\D) is of even par- 
ity, and the triplet t]^ n\-^) i^ ^^ ^'^^ parity. Apart from 
parity the two-hole states can be classified additionally 
according to their total spin, with four subspaces, i.e., 
5 = even, S — odd, S — 1 even, and S — 1 odd aris- 
ing. Two-pseudo-fcrmion scattering induced by Hamil- 
tonian ( [l8[pi| ) occurs only in the the 5 = even channel 
(via Hs) and in the 5 = 1 odd channel (via Ht)- There- 
fore we will focus on these channels in the following. 



A. 5* = even channel 

To calculate the even parity 5 = energy spectrum 
we start from the hole-pair Green's function 

G,{k,t) ^ ~zQ{t){D\[b^{t)bl]\D), (27) 

where the bare hole-pair propagator is set by Hi, of ( pi| ) 
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FIG. 2. T-matrix hole-pair Green's function Gb{k,uj). 
Squares: scattering vertex Skq (|l9|). Circles: hard-core re- 
pulsion Uq (psi). Full lines: SCBA pseudo-fermion Green's 
functions. Dotted lines: bare hole-pair Green's functions. 
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and for the remainder of this work retarded propagators 
with w = a; + 477 and 77 — > 0+ are implied. Following the 
preceding discussion the interacting propagator Gb{k,uj) 
is calculated by the T-matrix summation of fig. 



Gb{k,uo) = {Q Gl{uo)){U-M{k,uo)y 
where M{k, w) is a 3 x 3 matrix with 

Mii(22){k,U}) = — ^/oo(7r7r)(fc,'?,w) 
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with Skq from ( |28| ) and Uq from (25). Ifi,^{k,q,uj) refers 
to the two-pseudo-fermion bubble 

I^u{k,q,io) = / duj'{lmGf,^a{q-,^')) 



G,Ak-q,^-^') (31) 

with /i = 0, TT, u = 0, TT and the SCBA renormalized 
pseudo-fermion Green's functions Gky,a{k, uj) of (vW- On 
the l.h.s. of ( pl|) the spin index a has been omitted since 
Gky.<j{k,ijj) and thus I^^(k,uj) is spin-independent. This 
implies that, while collective excitations may appear only 



in the S — even and the S = 1 odd subspace, the two- 
pseudo-fermion continua are degenerate with respect to 
total spin. For even parity states only diagonal elements 
of I^i, enter the calculation, odd parity states imply off- 
diagonal elements. 



B. 5 = 1 odd channel 

The odd parity part of the spectrum is extracted from 
the triplet-hole-pair Green's function 
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where 72, , is the Bogoliubov-transform (|5|) of the triplet 
and the choice of the z triplet is arbitrary due to 
spin rotational invariance. The form factor g^q — 
{2/NY/'^UqC0s{^ — q) has been introduced for com- 
putational convenience: the matrix elements of Ht 
which enter the T-matrix resummation are of the form 
UqUq'Tkqq' = '2,^/2t2gkq9kq' whcrc, duB to the limit Ji ^ 
J2, we have neglected contributions oc Vq. Thus guq re- 
duces the number of distinct two-particle reducible el- 
ement to the T-matrix summation. This form factor 
changes the weight but not the position of the T-matrix 
eigenvalues. 

The hole-pair-triplet Green's function is coupled to 
two pseudo-fermion scattering states via Ht- With re- 
spect to the interactions Ht and U the T-matrix evalua- 
tion, fig. 3, proceeds analogous to that in 5* = even 
subspace. Note however, that to the zeroth order in 
the pseudo-fermion scattering states the hole-pair-triplet 
Green's function G1^{k,u}) 

GUk,uj) = -^Y.9lqJdco'{lmD{q,Lo') x 

Gb{k-q,uj- uj') = ^ 9lqGb{k - q,uj - lUq) (33) 



FIG. 3. T-matrix hole-pair-triplet Green's function 



Gbt(k,uj). Squares: Scattering vertex Tj 



kqq' 



Circles: 



hard-core repulsion C/o (E5|). Dashed line: Triplet Green's 
function. Full line: SCBA pseudo-fermion Green's function. 
Dotted line with box: T-matrix hole-pair Green's function of 
fig. I 




FIG. 4. Hole-pair-triplet spectral function Aijt{k,UL>) at 
fc = 0. Ji/fi = 3, J2/ti = 0.3, tx=t2 = l, N = 32. 
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(35) 



contains interaction effects already, since Gb{k,uj) is the 
interacting hole-pair Green's function ( p9| ) of Fig. y. 
D{q,uj) = l/(w — LOq) is the triplet Green's function ob- 
tained from (0). The full T-matrix reads 



GM{k,u;)^{0 G',,{k,u;)){h~ M'{k,u;jr' f^^') 



with the 2x2 matrix 
U, 



^ii = l^E^o.(A;,g,^) 
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Gbtik,^)^gkqIo7v{k,q, 



(34) 
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V. RESULTS AND DISCUSSION 

The energy spectrum in the even (odd) parity sector 
of the two-hole spectrum results from an evaluation of 
the hole-pair- (triplet) Green's functions Gb(bt){k,uj) and 
a subsequent determination of the regions of nonvanish- 
ing spectral density Ab(j,t)(k,Lo) = —ln\Gb(bt){k,uj)/T:, an 
example of which is depicted in Fig. 0^ for fc — at 
Ji/ti = 3, J2/ii = 0.3, ti = t2 = 1, and N = 32. In 
evaluating Gb(bt), both the lo integrations and k sum- 
mations are performed numerically for small, but finite 
77 and Uo/ti — 10^. We have checked that beyond the 
latter value the spectrum is insensitive to any further 
increase of C/q- Moreover, and except for later compar- 
ison to results of exact diagonalization we use a system 
size of A^ = 32 rungs for the remainder of this work. We 
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FIG. 5. Low lying two-hole excitations of a spin ladder with 
A'' = 32 rungs and Ji/ti = 3, J2/ti — 0.3 and fi — t2. (a) 
even parity (b) odd parity. 



have checked that for larger system sizes the spectrum re- 
mains almost unchanged suggesting that our results are 
sufhciently close to the thermodynamic limit. 

Figure Ha (b) shows the dispersion of the even (odd) 
parity low energy two-hole excitations of a spin ladder 
with J2/ti — 3, Ji/ii = 0.3 and ^2 = ti- The state of 
lowest energy is a two-hole bound state with 5 = and 
even parity, which is split off from the continuum. The 
binding effect is due to two holes residing on the same 
rung and thereby maximizing the number of singletstJ. 
The overall ground state of the ladder is this S = even 
bound-state at fc = 0. The charge excitations of low- 
est energy are intra-band excitations in the S = Q even 
bound state and therefore gapless, while the lowest spin 
excitations are inter-band excitations between the ground 
state and the S = 1 even coptinuum and therefore are 
gapful— .Thus a Luther-EmeryCHI rather than a Luttinger- 
liquidE^ applies to the present case. This ground state is 
consistent with that found by Lanczos diagonalizatjonfl'El, 
density matrix renormalization group (DMRG)I13, and 
high order series expansioncZI. 

In the odd parity spectrum of fig. pp an S* = 1 bound 
state is observed which- is-|Consistent with results from 
exact diagonalizationErEil~E3. In the T-matrix this state 
arises as a collective excitation due to local mixing of 
the hole-pair-triplet state with two pseudo-fermions. For 
very strongj-epupling Ji ^ J2,ti, ^2 an additional anti- 
bound stateE^ arises, which however merges with the con- 
tinuum already at intermediate couplings, because the 
continuum is strongly spread due to the incoherent part 
of the single-hole spectra. 

To emphasize the difference between the bound states 
in the 5* = and S = \ sector, we show the binding en- 
ergy Eb = min(£^cont) — niin(ii^bs) as a function of Ji for a 




FIG. 6. Bound state binding energy Eb for 5* = (squares) 
and 5* = 1 (diamonds). J2/J1 = 0.1, fi = t2, N = 32. Inset: 
maximum k at which the bound state exists. 



constant ratio J2/J1 =0.1 in fig. |^ where min(i?(.ont(bs)) 
denotes the minimum energy of the continuum (bound 
state). For large Ji a linear behavior of the S = binding 
energy is found in contrast to a constant binding energy 
for the S — I state. This difference reflects the differ- 
ent binding mechanisms: the energy gain of the S = 
bound state is due to. the additional singlet which ac- 
companies this stateEj. In this case the energy gain is 
ex Ji. In the S = 1 sector, and regarding Ji only, the 
local hole-pair-triplet state is energetically equivalent to 
a two pseudo-fermion scattering state. Therefore bind- 
ing in this channel is mainly due to kinetic delocalization 
of the hole-pair-tHplet into two pseudo-fermions on the 
neighboring sitesEJ. We note, that only for a finite range 
of parameters Ji/ii ^ 1 the bound states exist over the 
entire Brillouin zone. This is shown in the inset of fig. 
|6| which depicts the maximum values of k for which the 
bound states can be separated from the continuum. 

Returning to fig. ||(a) and (b) it is evident, that in 
the even parity sector the high energy states form two 
continua, degenerate in 5 = and S — I. The lower 
one of these, i.e. between E « — 4ii and E « — ti, is 
due to two pseudo-fermion scattering states formed out 
of the bonding orbitals, while the higher ones are due to 
either the incoherent parts of the bonding spectra or to 
two pseudo-fermion scattering states both in antibonding 
orbitals. Note that in fig. || we have cut oS^ the spectrum 
at high energies with additional continua existing above 
the cut off. 

The lower bound of the S = 1 odd continuum in 
fig|q(b) shows two local minima at fc = and k = tt. 
This structure is due to the twp-.types of excitations 
which contribute to the continuuirO: (i) scattering states 
of a triplet and the hole-pair on a rung, i.e. the hole- 
pair-triplet bubble in fig. 3 of the diagrammatic theory, 
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FIG. 7. Low lying two-hole excitations of a spin ladder 
with N ^ 8 and Ji/ii = 3, Jg/ti = 0.3 and ti = t2. Thin 
lines: reproduced after Troyer et al.^. The perturbative spec- 
trum has been shifted by a constant ground state energy of 
Eo = -ig.lii. 



and (ii) the continuum of two pseudo-fermions, i.e. the 
pseudo-fermion bubble in fig. 3. The latter is degenerate 
with the S — odd excitations and its minimum of en- 
ergy is at fc = TT, while the former one has no counterpart 
with S = and has its minimum at fc = 0. This excita- 
tion has a dispersion of its lower edge which is reminis- 
cent of that of the single triplet on the undoped ladder. 
For the parameters chosen in Fig. (ph both minima of 
the continua are nearly degenerate whereas in general it 
depends on the particular choice of ti,t2, Ji, J2, which 
minimum is at the lowest energy. 

Next a comparison of our diagrammatic results with 
exact diagonalization data is performed. To this end, in 
fig. ^ a twoJiole spectrum reproduced from the work of 
Troyer et al.B is depicted together with our results from a 
T-matrix evaluation. In order to perform the comparison 
identical parameters, and in particular identical systems 
sizes have been chosen. Thick lines in this figure repre- 
sent the diagrammatic, thin lines the numerical result. 
Obviously the agreement is very good. We emphasize 
that this kind of agreement is promoted by the limit of 
strong rung coupling J1/J2 ^ 1 which pertains to fig. 0. 
We do not expect similar agreement for J2I Ji -^ 1- 

To conclude this section we discuss the spin-gap of the 
doped ladder. From the spectrum of fig. ^ it is obvious 
that various S = 1 states are found in the doped ladderO: 
(i) the even parity continuum, (ii) the odd parity bound 
state, and (iii) the odd parity continuum consisting (a) 
of triplet-hole-pair scattering states with a minimum at 
k = -K and (b) of two pseudo-fermion scattering states 
with minimum at A: = 0. All of these excitations are 
gapped, with the actual spin gap set by the lowest energy 
S = 1 excitation. Fig. || shows the excitation energies for 
all of the different S = 1 states for Ji/ti = 3 and ii = ^2 
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FIG. 8. Excitation energies of the different 5 = 1 states as 
function of J2. Ji/fi = 3.0, h = t2, N ^ 32. 



as a function of J2 . The spin gap of the undoped ladder is 
also displayed for comparison. For nearly the entire range 
of J2 in fig. H the even parity continuum (i) defines the 
spin gap. The gap-size is set by the energy difference to 
the 5 = bound state at fc = 0. The continuum (i) has 
no counterpart in the undoped ladder. Thus we find that 
the spin gap evolves discontinuously upon doping in this 
parameter regime. The excitations at higher energies are 
the bound state (ii) and the continuum (iii) for J2/ti < 
0.9, while the bound states merges with the continuum 
at this particular value of J2. For J2 > 0.3 the lower 
edge of the odd parity continuum is set by the triplet- 
hole-pair scattering state, i.e. (iii a). As noted above 
this edge of the continuum has a dispersion similar to 
that of the triplet excitation on the undoped ladder and 
therefore leads to a spin gap also comparable to that of 
the undoped ladder. For J2/ti < 0.3 the lowest odd 
parity continuum state is the two pseudo-fermion state, 
i.e. (iii b), which explains the deviation from the spin 
gap of the undoped ladder. Finally at J2/ti « 1.2 a level 
crossing occurs and the hole-pair-triplet continuum, i.e. 
(iii b), becomes the lowest 5=1 excitation. Only in this 
parameter regime the spin gap evolves continuously upon 
doping. 

Unfortunately the LHP approximation considerably 
underestimates the spin-gap of the triplet excitation for 
J2 — > Ji compared to numerically exact resultaSJ. There- 
fore the results for the spin-gap in Fig. g are qualitative 
only as J2 increases. Nevertheless the level crossing cited 
above has been found also in numerical diagonalizationE3, 
however at a parameter range not accessible to our dia- 
grammatic approach. Introducing a diagonal next-near- 
neighbor hopping t' between site 1(2) of rung n and site 
2(1) of rung n+1 it has been shown in numerical studiesc2l 
that the spin gap can be tuned continuously from being 
set by the triplet of the undoped ladder via the holc-pair- 
triplet bound-state to the even parity continuum. Fi- 
nally, very recently, it was found that a strongly negative 



t' changes the nature of the ground state from the-|5' = 
bound state to that of the SL^ even continuurrO, thus 
from a Luther-Emery hquidE3 to a scenario with gapless 
charge and spin excitations. 



VI. CONCLUSIONS 

In conclusion we have analyzed the excitations in spin 
ladders doped with two holes. Applying a mapping of 
the generalized t-J model onto a coupled boson-fermion 
model the impact of both, scattering by the spin back- 
ground and hole-hole interactions on the two hole spec- 
trum has been analyzed in detail. We find the ground 
state to be a two-hole bound state, whereas excitations 
consist of various continuum states with S = 0,1 and 
even or odd parity. Additionally, depending on the par- 
ticular choice of parameter, S — 1 odd parity bound and 
antibound states are found. The binding energy of the 
bound states and the spin gap of the system have been 
analyzed. Our results compare well with numerical stud- 
ies of finite systems. Yet, our method allows for a study 
of system sizes close to the thermodynamic limit. 
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APPENDIX A: 

In this appendix, and for the sake of completeness we 
list the Hamiltonian Hi of (p^, i.e. of the single- hole 
sector 



Hi = -h Y^ aL„,^a2^„ ,, + h.c. 



+ t2S 



5]4(-iri(sri,„ + s^^i,„) + /i., 






5]ti(-iri s^-Li,,_i + s4^i^„+i +/. 



],ri 



- i2 XI * (4 X t„+i) S^/+i „ + h.c. 
j-n 

■ 2_^ i (t„ X t„) (^S„^;^ „_|_j + S„_;^ „_;^ j . (Al) 



2 



where tj^ = (ij; „, il „, t^ „). This Hamiltonian includes 
(i) intra-rung hopping, (ii) inter-dimer hopping and (iii) 
exchange scattering. Interdimer hopping can be ei- 
ther spin-diagonal or accompanied by spin-flip scatter- 
ing. This includes (ii,a) singlet-singlet and (ii,b) singlet- 
triplet and (ii,c) triplet-triplet transitions of the spin 
background upon hole doping. For further details see 
ReiM. 



J,n 



1 Pof reviews see: E. Dagotto and T. M. Rice, Science 271, 
618 (1996); M. Takano, Physica C 263, 468 (1996); S. 
Maekawa, Science 273, 1515 (1996); T. M. Rice, Z. Phys. B 
103, 165 (1997); H. Tsunetsugu, Physic a B 237-238, 108 
(1997). E. Dagotto, |cond-mat/9908250| ; Rep. Prog. Phys. 
62, 1525 (1999). 

^ Z. Hiroi, M. Azuma, M. Takano and Y. Bando, J. Solid 
State Chem. 95, 230 (1991); M. Azuma, Z. Hiroi, M. 
Takano, K. Ishida and Y. Kitaoka, Phys. Rev. Lett. 73, 
3463 (1994). 

^ H. Iwase, M. Isobe, Y. Ueda and H. Yasuoka, J. Phys. Soc. 
Jpn. 65 2397 (1996). 

* M. Uehara, T. Nagata, J. Akimitsu, H. Takahashi, N. Mori, 
and K. Kinoshita, J. Phys. Soc. Jpn. 65, 2764 (1996). 

^ E.M. McCarron, M.A. Subramanian, J.C. Calabrese, und 
R.L. Harlow, Mater. Res. Bui. 23, 1355 (1988) 

® H. Tsunetsugu, M. Troyer, T. M. Rice, Phys. Rev. B 49, 
16078 (1994). 

■^ H. Tsunetsugu, M. Troyer, and T. M. Rice, Phys. Rev. B 
51, 16456 (1995). 

* M. Troyer, H. Tsunetsugu, T. M. Rice, Phys. Rev. B 53, 
251 (1996). 

^ S. Haas and E. Dagotto, Phys. Rev. B 54, R3718 (1996). 
^° S. R. White, D. J. Scalapino, Phys. Rev. B 55, 6504 (1997). 
^^ M. Brunner, S. Capponi, F. F. Assaad, A. Muramatsu, 

Phys. Rev. B 63, 180511(R) (2001). 
^2 R. Eder, Phys. Rev. B 57, 12832 (1998). 
" O. P. Sushkov, Phys. Rev. B 60, 3289 (1999). 
" C. Jurecka, W. Brenig, Phys. Rev. B 63, 094409 (2001). 
^^ E. Dagotto, J. Riera, D. Scalapino, Phys. Rev. B 45, 5744 

(1992). 
^^ C. Gazza, G. B. Martins, J. Riera, E. Dagotto, Phys. Rev. 

B 59, 709 (1999). 
" T. Siller, M. Troyer, T. M. Rice, S. R. White, Phys. Rev. 

B 63, 195106 (2001). 
^** T. F. A. Miiller, T. M. Rice, Phys. Rev. B 58, 3425 (1998). 
" S. Rommer, S. R. White, D. J. Scalapino, Phys. Rev. B 61, 

13424 (2000). 
^° D. Poilblanc, D. J. Scalapino, W. Hanke, Phys. Rev. B 52, 

6796 (1995). 
^^ E. Dagotto, G. B. Martins, J. Riera, A. L. Malvezzi, C. 

Gazza, Phys. Rev. B 58, 12063 (1998). 
^^ J. Riera, D. Poilblanc, E. Dagotto, Eur. Phys. J. B 7, 53 

(1999). 
^^ D. Poilblanc, O. Chiappa, J. Riera, S. R. White, D. J. 

Scalapino, Phys. Rev. B 62, R14633 (2000). 



K. Tsutsui, D. Poilblanc, S. Capponi, |cond- mat/0106389 



C. A. Hayward, D. Poilblanc, R. M. Noack, D. J. Scalapino, 
W. Hanke, Phys. Rev. Lett. 75, 926 (1995). 

^^C. A. Hayward, D. Poilblanc, Phys. Rev. B 53, 11721 

(1996). 
^'^ J. Oitmaa, C. J. Hamer, Z. Weihong, Phys. Rev. B 60, 

16364 (1999). 
^** M. Sigrist, T. M. Rice, F. C. Zhang, Phys. Rev. B 49, 

12058 (1994). 
^^ Y. L. Lee, Y. W. Lee, C.-Y. Mou, Z. Y. Weng, Phys. Rev. 

B 60, 13418 (1999). 
^° A. Lauchli, Diploma thesis, ETH Ziirich (1999). 
^^ H. Lin, L. Balents, M. P. A. Fisher, Phys. Rev. B 58, 1794 

(1998). 
^^ I. Bose and S. Gayen, J. Phys.: Cond. Mat. 11, 6427 

(1999). 
^^ G. Sierra, M. A. Marti'n-Delgado, J. Dukelsky, S. R. White, 

D. J. Scalapino, Phys. Rev. B 57, 11666 (1998). 

^* S. Sachdev and R. N. Bhatt, Phys. Rev. B 41, 9323 (1990). 
^^ A. V. Chubukov, Pis'ma Zh. Eksp. Teor. bf 49, 108 (1989); 

[JETP Lett. 49, 129 (1989)]. 
^'^ A. V. Chubukov and Th. Jolicoeur, Phys. Rev. B 44, 12050 

(1991). 
^^ O. A. Starykh, M. E. Zhitomirsky, D. I. Khomskii, R. R. 

P. Singh, and K. Ueda, Phys. Rev. Lett. 77, 2558 (1996). 
^* W. Brenig, Phys. Rev. B 56, 14441 (1997). 
■^^ T. Barnes, E. Dagotto, J. Riera, and E. S. Swanson, Phys. 

Rev. B, 3196 (1993). 
■*'' M. Reigrotzki, H. Tsunetsugu, and T. M. Rice, J. Phys.: 

Condens. Matter 6, 9235 (1994). 
*i C. Jurecka and W. Brenig Phys. Rev. B 61, 14307 (2000). 
*^ Very recently Park and Sachdev used this representation 

for a mean-field theory of two-dimensional antiferromag- 

nets and spin ladders at finite doping. 

K. Park, S. Sachdev, |cond-mat/010451£| (2001). 



43 



O.P. Sushkov and V.N. Kotov Phys. Rev. Lett. 81, 1941 

(1998). 
■** V. N. Kotov, O. P. Sushkov, Z. Weihong, and J. Oitmaa, 

Phys. Rev. Lett. 80, 5790 (1998). 
''^ V. N. Kotov, O. P. Sushkov, and R. Eder, Phys. Rev. B, 

6266 (1999). 
''^ C. Jurecka, V. Griitzun, A. Friedrich, W. Brenig, Eur. 

Phys. J 21, 469 (2001). 
*^ S. Schmitt-Rink and CM. Varma, Phys. Rev. Lett. 60, 

2793 (1988). 
** C. L. Kane, P. A. Lee, N. Read, Phys. Rev. B 39, 6880 

(1989). 
"^ G. Martinez, P. Horsch, Phys. Rev. B 44, 317 (1991). 
^° Z. Liu, E. Manousakis, Phys. Rev. B 45, 2425 (1992). 
^^ S. Gopalan, T. M. Rice, M. Sigrist, Phys. Rev. B 49, 8901 

(1994). 
^^ A. Luther and V.J. Emery, Phys. Rev. Lett. 33, 589 (1974); 

V.J. Emery, in Highly Conducting One-Dimensional Solids, 

edited by J.T. Devreese et al. (Plenum, New York, 1979). 
^^ F.D.M. Haldane, Phys. Rev. Lett. 45, 1358 (1980); J. Phys. 

C 14, 2585 (1981). 
^* M. Greven, R. J. Birgeneau, and U. -J. Wiese, Phys. Rev. 

Lett. 77, 1865 (1996). 



